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ABSTRACT 

X-ray Flash (XRF) 100316D, a nearby super-long under-luminous burst with a peak energy Ep ~ 20 
keV, was detected by Swift and was found to be associated with an energetic supernova SN 2010bh. 
Both the spectral and the temporal behavior are rather similar to XRF 060218, except that the latter 
was associated with a "less energetic" SN 2006aj, and had a prominent soft thermal emission compo- 
nent in the spectrum. We analyze the spectral and temporal properties of this burst, and interpret 
the prompt gamma-ray emission and the early X-ray plateau emission as synchrotron emission from 
a dissipating Foynting flux dominated outflow, probably powered by a magnetar with a spin period 
of P ~ 10 ms and the polar cap magnetic field i?p ^ 3 x 10^^ G. The energetic supernova SN 2010bh 
associated with this burst is however difficult to interpret within the slow magnetar model, and we 
suspect that the nascent magnetar may spin much faster with an initial rotation period ^ 1 ms. It 
suggests a delay between the core collapse and the emergence of the relativistic magnetar wind from 
the star. The diverse behaviors of low-luminosity GRBs and their associated SNe may be understood 
within a unified picture that invokes different initial powers of the central engine and different delay 
times between the core collapse and the emergence of the relativistic jet from the star. 
Subject headings: Gamma rays: general — Radiation mechanisms: non-thermal — X-rays: general 



1. INTRODUCTION 

After four years of waiting since the detec- 
tion of X-Ray Flash (XRF) 060218/SN 2006aj, an- 
other pair of low-luminosity (LL) GRB - super- 
nova (SN) association, XRF 100316D / SN 2 010bh 
at redshift z = 059 (IVergani et al. I 120100 . was 



captured by 
16, 2010 



Swift (IGehrels et al. 
(iStamatikos et al. I 12010: 



12001 



on March 



Wiersema et al. 



\2Um iChornock et al. I 120101: iRau et al. I I20100 . with 
a detection rate fully consistent with the pop- 
ulation studies of these nearby LL - GRB events 
(ICoward I 120051: ISoderberg et al. I 120061: ILiang et al.1 
I2007at IGuetta fc Delia Valle 1 120071) . Before this event, 
four pairs of nearby {z < 0.2) secure GRB(XRF)- 
SN associations have been identified. These are GRB 
98042 5/SN 1998bw at z = 0.0085 (e.g., IGalama et ahl 
1991. GRB 03032 9 /SN 2003dh at z = 0.168 (e.g., 
Hiorth et al.1 [200l. GRB 031203/SN 20031w at z 



0.105 l'e.g.. lMalesani et al"~|[2004). and XRF 0602 18/SN 
2006aj at z = 0.0331 (e.g.. ICampana et al. |[2006[) . The 



nature of the GRB/SN connection and the interplay be- 
tween the GRB and the SN components are still poorly 
understood. In this paper, we analyze and interpret the 
Swift BAT and XRT data of XRF 100316D, paying spe- 
cial attention to the similarities and differences between 
the XRF 100316D/SN 2010bh and XRF 060218/SN 
2006aj. 

2. THE 7-RAY AND X-RAY OBSERVATIONS 

GRB 100316D triggered Swift/BAT at 12:44:50 UT 
(Ttrig), and the XRT began observing the field 137.7 sec- 
onds after the BAT trigger (Stamatikos et al. 2010). A 

|yzfan®pmo. ac.cn (YZF), zhang@physics.unlv.edu(B"Z)] 



bright, steady un-catalogued X-ray source was detected. 
By analyzing the BAT survey data before and after the 
trigger time, it is found t hat this event possibly sta rted at 
~ (Ttrig - 1500) seconds ([Stamatikos et al. Il2010l) . How- 
ever, the gamma-ray fiux of the source kept almost con- 
stant from Ttrig — 1500 to Ttiig — 500, before increasing 
significantly at Ttrig — 500. We therefore take the starting 
time of this event as Tti-is — 5 00 seconds. Th e host galaxy 
redshift is z = 0.059 (.Vergani et al. |[2010[ ). An associ- 
ated supernova, SN 20i0bh, was detected by Wiersema et 
al. (2010) and spectroscopically confirmed by Chornock 
et al. (2010). 

We extract the observed lightcurves and spectra 
of XRF 100316D from the BAT and XRT event data. 
The details of ou r da ta reduction are pr esented in 
IZhang et all (|2007D and lLiang et al. I ()2007bn . The joint 
BAT/XRT spectrum from Ttng + 138 s to Ttrig + 736 s 
is well fit with a Cutoff Power-law model, with a power- 
law index of 1.32 ± 0.03 and a spectral peak energy Ep 

of I9.6tli keV (xVdof = 1320/1057), Fig^. This 
suggests that the X-ray and the gamma-ray emissions 
are from the same emission component. A time-resolved 
spectral analysis during the period from Ttrig + 138 s to 
Ttrig + 736 s shows that Ep clearly evolves with time, 
from 32.6lg;'5^ keV (from Ttrig + 138 s to Ttrig + 240 s) 

to IS.Stli keV (from Ttrig + 240 s to Ttrig + 734 s). The 
un-absorbed lightcurves and the temporal evolution of 
Ep are shown in Figure [H along with the data of GRB 
060218. The late X-ray lightcurves of the other nearby 

^ The time-integrated BAT spectrum can be reasonably fitted by 
a single power <x i/-2-5±0.3 Together with the joint BAT-XRT 
fit, we are able to constrain Ep ~ 20 keV and therefore define this 
burst as an XRF. 
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GRBs, including GRB 980425, 031203, and 030929, are 
also displayed in Figure [2j The X-ray light curves of 
XRF 100316D and XRF 060218 are rather similar. Both 
have an plateau phase extending to more than 1000 s. 
In the case of XRF 100316D, there is an observational 
gap between Ttrig + 736 s and Ttrig + 36437 s. The late- 
time observational data (after Ttrig + 36437 s) is soft and 
consistent with a single power law decay with a decay 
index a = 1.31±0.21. This is similar to the late X-ray 
lightcurve of XRF 060218. XRF 060218 also shows a 
steep decay phase between the plateau and the late sin- 
gle power law decay phase. Although this is not observed 
in XRF 100316D, the data before Ttdg + 734 s and after 
Ttrig + 36437 s are consistent with having such a steep 
decay segment in between. 

There are however some differences between the two 
events. First, the X-ray data of XRF 060218 demand a 
soft thermal emission component with an evolving tem- 
perature, which has been interpreted as due to shock 
break out from the progenitor star (jCampana et al.~l 
120061 ). For XRF 100316D, we find that the data 
do not demand such a component. We noticed that 
iStarling et al. I ()2010( ) claimed a thermal component with 
kT ^ 0.14 keV and an energy that is ^ 3% of the entire 
X-ray emission. To check the consistency, we fit the time 
dependent BAT-I-XRT joint spectra by a Cutoff Power- 
law model (with absorption from both Milk Way and 
the host galaxy, wabs*zwabs*cutoffpl in Xspec 12) and a 
Blackbody-|-Cutoff Power-law model. We find that both 
models can give equally acceptable fits to the data. For 
example, for the slice 2 (from Ttrig -I- 240 s to Ttrig -I- 734 s), 
the former model gives jdof — 1305/1085 with A^gai = 
7.05 X 10^° cm-2 and iVH,iiost = 4.6 x lO^^cm^^ while 
the latter gives /dof ~ 1264/1085 with a larger host 
galaxy absorption iVH,host — 1-3 x lO^^ci n"^ and a ther- 
mal te mperature that is consistent with IStarling et al. I 
(|201Cf) . Note that the peak of the proposed black body 
component is near the low end of the XRT energy band 
and the neutral hydrogen absorption around the thermal 
peak is large. This makes great uncertainty on iden- 
tification of such a thermal emission component from 
the data. Our spectral analysis cannot robustly identify 
this component, although the is slightly improved by 
adding it. Therefore, we do not claim a thermal emission 
in the observed spectrum and only stick to the Cutoff 
power-law model to discuss possible theoretical implica- 
tions. The second difference between the two events lies 
in the supernova data. The modeling of SN 2006a.i sug- 
gests a kinetic energy ^ 2.5 x 10^^ erg fM azzali et al. I 
[2006b), much smaller than that of other SNe associated 
with nearby GRBs (see Fig 13]). Although the modeling 
of SN 2010bh is not available yet, curre nt data imply an 
SN e vent as energetic as SN 1998bw (jChornock et al. I 
l2010f R. which is about one order of magnitude more en- 

^ With the simplest assumptions that (1) the opacity of the 
SN outflow is from Thompson scattering of electrons and is 
~ 0.2 cm^ and (2) The density of expanding SN material 
takes the form oc R^^ (fc ~ 6 — 8, J. S. Deng, private com- 
munication, s ee also Matzncr & McKcc 1999; Bcrgcr et al. 2002; 
IChevalier fc Fransson .■200S ^ , it is straightforward to show that the 
mass of the SN material moving faster than Va is Mg^{> Vs) ~ 
Stt^^ Vg^t^, where Vb is the photospheric velocity at t, a^Y is 
the Thompson scattering cross section and nip is the rest mass 






Fig. 1.— The joint BAT-XRT spectrum of XRF 100316D with a 
Cutoff Power-law model fit (line) for the time-integrated spectrum 



from Ttr 



138 s to Ttr 



736 s. The model parameters are 



1.32 ± 0.03 and En 



19.6 



+3.3 



keV. 
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Fig. 2. — Upper: Unabsorbed luminosity lightcurve in the XRT 
band of GRB 100316D in comparison with GRBs 980425, 030329, 
031203, and 060128. The BAT data of GRB 100316D are extrap- 
olated to the XRT band. Lower: Comparison of the Ep temporal 
evolution of GRB0100316D with GRB 060218. The Ep data of 
GRB 060218 are taken from Toma et al. (2007) 

ergetic than SN 2006aj. So from the observational point 
of view, XRF 100316D/SN 2010bh and XRF 060218/SN 
2006aj are not strictly "twins" . 

of proton. The corresponding kinetic energy is -Ek,SN(> Ki) ^ 



' k^^K*t'^- For SN 2010bh, Vs ~ 0.09c at t ~ 21 day, we have 
£^k,SNiObl(> 0.09c) ~ 1.6 X 1052 [(k - l)/3(fc - 5)] erg. 
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Fig. 3. — The isotropic energy of the prompt emission vs. the 
kinetic energy of the supernova outflow. The kinetic energy of SN 
2010bh is estimated to be lar ger than ~ 10^^ erg (see footnote 
\2f. Other data are taken fromju] I^OOB ). The possible maximum 
energy ~ 5 X 10^^ erg that can be provided by a pulsar with P < 1 
ms and / ~ 2 X 10^^ g cm^ is also plotted. 

3. A POSSIBLE MODEL FOR THE LONG-LASTING X-RAY 
PLATEAU 

In the following, we define t = T — Ttrig + 500 s, 
i.e. the time elapse since Ttri<y — 500 s. We interpret 
all the BAT/XRT data of XRF 100316D for < i < 
1.23 X 10^ s as "prompt emission" (i.e. the radiation 
powered by some internal energy dissipation processes) 
for the following two reasons. First, the steady plateau 
behavior observed in both BAT and XRT band at i < 
1.23 X lO'^ s with an evolving Ep is difficult to interpret 
within afterglow models. Second, the sharp decline of the 
X-ray emission or even steeper) expected in the time 
interval 1.23 x lO'^ s < i < 3 x 10'' s resembles the early 
rapid decline that has been detected in a considerable 
fraction of Swift GRBs, which is widely ta ken as a piece 
of ev idence of the end of prompt emission (jZhang et al. I 
[200^ . The nature of the X -ray emission detected at 
t > 3 X 10^ s is hard to pin down. Its spectrum is very 
soft (photon index T = 3.3^g'^), similar to that of XRF 
060218. This is also unexpected in the external forward 
shock models, and this late X-ray component may be 
related to a l ate central engine afterglow, whose origin is 
unclear fe.g.. iFan et"ani2006HSoderberg et al. II2006D . 

The prompt BAT/XRT data do not show a signifi- 
cant variability (Figl2]). The time-averaged 7— ray lumi- 
nosity is ~ 3 x 10^^ erg s~^ and the X-ray luminosity 
is ~ 2 X 10^^ erg s~^. The bolometric luminosity of the 
XRF outflow is therefore expected to be in the order of 
10^^ erg s~^. The duration of the BAT emission is at 
least 1.23 x 10'^ s, and can be longer. The relatively 
steady energy output is naturally produced if the central 
engine is a neutron star with significant dipole radiation. 
The dipole radiation luminosity of a magnetized neutron 
star can be described as 

Ldip = 2.6 X 10^^ erg Bl^^Rle^l (^1 + ,(1) 

where Bp is the dipole magnetic field strength of the 
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Fig. 4.— Broadband SED from UVOT, XRT, and BAT data. 
Grey points show the time-averaged BAT -f XRT spectrum between 
TtrigH-lSO sand Ttrig-l-744 s. The thick dashed line represents 
an absorbed broken power-law fit (i.e., wabs*zwabs*bknpower in 
XSPEC although the absorption components are not plotted here) 
to the BAT-I-XRT data leading to Fi = 1.42, F2 = 2.48 and 
^'broak = 16 keV. The solid line shows the same fitting as above for 
higher energy band but with an additional break at 1 keV, below 
which photon index is set to —2/3. UVOT observations are taken 
from Starling et al. (2010). The extinction in each filter has been 
corrected by adopting Emw{B — V) = 0.12 from the Milky Way 
and i?host(^ — V) = 0.1 from the host galaxy (Starling et al.2010; 
Chornock et al. 2010) and a Milky Way-like extinction curve forall 
bands (Pei 1992). 

neutron star at the magnetic pole, Rg is the radius of 
the neutron star, fl is the angular frequency of radi- 
ation at i = 0, To = 1.6 X 10^B~^^^n:['^Ii5R~^ s is 
the corresponding spin-down timescale of the magne- 
tar, and / ^ 10**^ g cm^ is the typical moment of in- 
ertia of the magnetar (Pacini 1967; Gunn & Ostriker 
1969). Here the convention = Q/IO" is adopted in 
cgs units. One then has L^ip ^ const for t <^ tq and 
L(jip oc t^^ for i 3> Tg. An abrupt drop in the X-ray 
flux with a slope steeper than may be interpreted as 
a decrease of radiation efhciency, or the collapse of the 
neutron star into a black hole, possibly by losing the an- 
gular momentum or by accreting materials. Within such 
a model, the fact that 

idip - lO"^ erg s"\ To - 1000 s, 

would require {Bp^i4, Q4, iJg.e) ^ (30, 0.06, 1, 1). 
This is a slow (P ~ 10 ms) magnetar (Bp ~ 3 x lO^^ G). 

The composition of this spindown-powered outflow 
is likely Poynting-flux-dominated. Besides the magne- 
tar argument (which naturally gives a highly magne- 
tized outflow), another argument would be the lack of 
a bright thermal component with a temperature kT ^ 
10 keV l\j'^RqI^^ from the outflow photosphere as 
predicted in the baryonic outflow model, where Rq is 
the initial radius wh e re the outfl ow is accelerated (e.g. 
IZhang fc Pe'eF|[2009l: iFa^lMl . One may argue that 
the photosphere radiation peaks at the observed Ep. 
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Such a scenario, however, is hard to account for the X-ray 
spectrum F^, oc iy-0-32±o.Q3 _ 

Below we discuss a possible magnetic dissipation sce- 
nario that would interpret XRF 100316D prompt emis- 
sion. By comparing the pair density (cx r is the 
radial distance from the central source) and the den- 
sity required for co-rotation (oc beyond the light 
cylinder of the compact object), one can estimate the 
radius at which the MHD condition breaks down, which 

reads tmhd 5 x 10^^ Llij.^ ai lty^rn,-2^r,i.5 where a 
is the ratio of the magnetic energy flux to the particle 
energy flux, Fi is the bulk Lorentz factor of the outflow, 
tv,m -P is t he minimum variability time scale of the cen- 
tral engine dZhang k Meszaros 2002 : Fan et al. 1 120051: 
iGao &: Fan 1 120061 ). Beyond this radius, intense elec- 
tromagnetic waves are generated and outflowing parti- 
cles a re accelerated (jUsov 11199^ iLvutikov fc Blackman I 
I2001f) . Part of the Foynting flux energy is converted to 
radiation. At tmhd , the comoving magnetic fields -Bmhd 
can be estimated as Bmhd ~ 20 ^cri,Qt~l^ _2 G, where 
^ < 1 reflects the efhciency of magnetic energy dissipa- 
tion. When magnetic dissipation occurs, a fraction Ec 
of the dissipated comoving magnetic energy would be 
eventually converted to the comoving kinetic energy of 
the electrons. Electrons may be linearly accelerated in 
the electric fields or stochastically accelerated in the ran- 
dom electromagnetic fields. One may assume that the 
accelerated electrons form a single power-law distribu- 
tion in energy, i.e. (in/d'ye oc 7"'' for 70 > 7e,mj where 
7o.m can be estimated as 70. m ^ 3.5 x 10^ (Ji.eCp, and 
Cp = (ee/0.5)(p— 2)/(p— 1). At rM Hn. the correspond ing 
synchrotron radiation frequency is (|Fan et al. 1120051 ) 



J^m,MHD ^ 2 X 10^^ Hz 'Ccri^ gCpFi^i.5t^^„ 



(2) 



The cooling Lorentz factor can be estimated by 
7c,c ~ 4.5 X fOi^Fi/(rMHDBMHD)> which gives 
a synchrotron cooling frequency i^c,mhd ^ 6 x 
10^'^ Hz Ff 1 5L4/^-3f7j:^tv,m,-2- The observed XRT 
spectrum can be approximated by F^^ oc j/-o.3~-o.4^ 
which is close to the fast cooling spectrum v^^^^, sug- 
gesting t'c.MHD S 10^^ Hz (it is straightforward to show 
that for the fiducial parameters adopted in this work the 
synchrotron self-absorption is below the optical band). 
On the other h and, the non-detectio n of prompt emis- 
sion by UVOT (jStarling et al. Il20f Of ) favors a high cool- 
ing frequency fc,MHD ~ 2.4 x fO^^ Hz. As shown in 
FigHl the synchrotron radiation model with such a high 
I'c.MHD can roughly account for the data. We then get a 
constraint 



Fi ^ 60Ll/'d(,'3^!.(fC,nl-2('^c,MHD/2.4 X fo" Hz)i 
Substituting this into eq.® we get 



/5 



54L 



-1/16 



i^c.MHD 



2.4 X IQi^ Hz 



wl/16c-l/2^-5/8.3/8 



Within such a scenario, Ep is defined by fm.MHD [Eq.([2|l]. 
An observed evolving Ep may be understood by invoking 
a decreasing tr, which is consistent with the magnetic 
field dissipation hypothesis. For i' > Ep/h (where h is 
Flanck's constant), the spectrum is oc v^^^"^ . 



The typical variability timescale of the radiation 
powered by the magnetic energy dissipation can be esti- 
mated as 

<Si - rMHD/(2r2c) ^ 10 s L^('(a/50)-ii,,„,_2(ri/60)-3, 

. . (3) 
which is much larger than the rotation period P 10 ms 
of the magnetar, implying that the pulsation of the cen- 
tral engine is undetectable. 

4. DISCUSSION AND SPECULATIONS 

So far the nearby supernova-associated GRBs, ex- 
cept GRB 030329, are found to be intrinsically under- 
luminous. They share the similarities such as low 
isotropic energies and smooth light curves, but differ 
in some aspects. For example, GRB 980425 and GRB 
031203 have shorter durations and higher Sp's than XRF 
060218 and XRF 100316D. The underlying physical pro- 
cesses that result in these differences are not well un- 
derstood. Through supernova modeling, it is found that 
GRB 980425 and GRB 031203 have a progenitor star 
massive enough to for m a black hole (jPeng et al. II2005I : 
iMazzah et aLll2006aD . while XRF 060218 has a less mas- 
sive proge nitor that most plausibly produces a neutron 
star (Maz zali et al. |[2006b[ ). The luminosity and the du- 
ration of XRF 100316D are consistent with the radia- 
tion from a neutron star with a dipole magnetic field 
-Bp ^ 3 X 10^^ G and a rotation period P ^ 10 ms. 
This seems to point towards a hypothesis that two types 
of central engines define the apparent dichotomy of the 
SN-associated LL-GRBs, i.e. black hole engines give rise 
to "shorter" and "harder" GRBs such as GRB 980425 
and GRB 031203, while magnetar engines give rise to 
very longand soft XRFs such as XRF 060218 and XRF 
100316DB 

Adding in SN data makes the scenario more compli- 
cated. Ahhough SN 2006aj associated with XRF 060218 
does not conflict with a slow magnetar central engine, 
SN 2010bh associated with XRF 100316D may be too 
energetic to be interpreted with a slow magnetar cen- 
tral engine. If it is confirmed that the kinetic energy 
of SN 2010bh is in excess of 10^^ erg, neither the neu- 
trino energy nor the magnetar spin energy ('^ lO'^'^ erg) 
are adequate to power the SN. A salient feature of the 
dipole spindown formula (Eq.[l]) is that if one shifts the 
time zero point (e.g. t' = t — to), the spindown law still 
applies, with the initial angular frequency re-defined as 
n' = ri(t' = 0), and the characteristic spindown time 

scale re-defined as Tq = 1.6 x W^E'll^W^^^ I45R^q. This 
suggests that the observed plateau feature can be still in- 
terpreted if the initial period is much shorter than 10 ms, 
say, Pq ~ 1 ms, if the time zero point is much earlier than 
Jtrig — 500 s (i.e. t = 0). This is because a power-law de- 
cay light curve may show an artificial plateau in the log- 
log space, if the zero time is mis-placed to a later epoch 
(|Yamazaki et al. II2009I: [Liang et al. Il2009[ ). Within such 
a scenario, a nascent magnetar was born with an initial 

^ We caution that such a scenario is not robust, since a magnetar 
engine may also drive SN 1998bw ( Wooslcy 2010), and a black hole 
engine may be also able to reproduce the XRF 100316D-like light 
curves through tuning the parameters of fall-back materials and ar- 
guing for a Poynting-flux-d ominated outflow from a highly magne- 
tized black hole engine (e.g. lMacFadven et al~l[200ll : [Zhang et al. I 
l200a) . 
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period Pq ^ 1 ms at t ^ — 5 x 10'^ s. Its intial dipole 
radiation was trapped by the envelope of the progeni- 
tor and could not escape. This spindown energy gives 
enough impetus to explode the star and power the ener- 
getic SN 2010bh. After a significant delay 5 x 10^ s 
to spin down from 1 ms to 10 ms for Bp ~ 3 x 10^^ G), 
the magnetar wind finally managed to escape as a rel- 
ativistic Poynting-flux-dominated outflow. An observer 
noticed the jet emission only around t — 0. The above 
argument also applies to the model of fallback accretion 
onto a nascent black hole. 

With such a hypothesis, one may envision a unified 
picture to understand the diversity of GRB/SN associ- 
ations, by invoking a variety of initial powers and the 
delay times between the core collapse and the emergence 
of the relativistic jet from the star. The speculation is 
the following: 

• To produce an energetic SN/luminous GRB (e.g. 
GRB 030329/SN 2003dh), the central engine is 
powerful (a black hole with an accretion disk or 
a rapidly spinning magnetar) and the relativistic 
outfiow can break out the progenitor soon enough 
when the engine is still working effectively. 

• To produce an energetic SN / underluminous GRB 
(e.g. GRB 980425/SN 1998bw, GRB 031203/SN 
20031w, and XRF 100316D/SN 2010bh), the cen- 
tral engine is initially powerful, but it takes time for 
the relativistic wind to emerge from the star. As it 
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breaks out the star, the central engine already fades 
down with a decreased power. The longer, softer 
XRFs are probably powered by a magnetar, while 
the shorter, harder GRBs are probably powered by 
a black hole. 

• To produce a less-energetic SN / underluminous 
GRB (e.g. XRF 060218/SN 2006aj), the central 
engine is a slow magnetar with an initial rotation 
energy less than 10^^ ergs. The emergence of the 
relativistic outfiow can be prompt or somewhat 
(but not significantly) delayed. 

Finally, a straightforward expectation from the spec- 
ulation that XRF 100316D outflow is Poynting-flux- 
dominated is that th e prompt emissio n should be lin- 
early polarized Ce.g.. iFan et al. |[2005[ ). The polarime- 
try measurements of events such as XRF 100316D and 
XRF 060218 would provide a criterion to differentiate 
this model from the shock breakout model, which does 
not predict a strong polarization signal. 
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